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SUMMARY 


This  contract,  transferred  to  the  auspices  of  the  SQUID  program 
effective  1  March  1981,  was  originally  planned  as  a  three  year  program. 

The  objectives  were  predominantly  experimental  in  nature: 
to  measure  the  unsteady  velocities  in  the  rotor  passage  of  a  single 
stage  axial  flow  compressor  and  map  the  velocity  field  at  operating 
conditions  near  rotor  blade  stall.  For  this  purpose  a  laser  doppler 
velocimeter  was  to  be  developed  and  married  to  a  new  single  stage 
axial  flow  compressor  designed  expressly  for  this  purpose. 

At  the  premature  conclusion  of  the  contract,  the  first 
successful  velocity  determination  in  the  running  compressor  had  not 
been  obtained,  although  the  LDA  and  its  associated  subsystem  was 
operational . 


Some  compressor  cascade  modelling  studies,  ultimately  intended 
to  be  guided,  or  "calibrated"  by  the  experimental  results  are  only/ 


partially  completed. 
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The  single  stage  compressor  was  constructed  under  this 


contract  and  a  variable  speed  control  procured  for  setting  the 
compressor  operating  point  in  conjunction  with  a  discharge  throttle 
valve.  This  system  is  running  satisfactorily  and  the  present  status 
may  be  described  as  being  in  the  midst  of  modifications  to:  a) improve 
the  optical  quality  of  the  laser  beam  penetration  and  return  through 
the  transparent  housing,  b)  mount  the  LDA  on  an  adjustable  frame  to 
permit  an  accurate  determination  and  viewer’  adjustment  of  the  measuring 


volume  location  in  3-space  and  c)  procure  and  install  of  a  set 
of  exit  guide  vanes  to  provide  more  realistic  compressor  loading 
characteristics  in  later  tests.  The  last  item  is  being  pursued  as 
part  of  a  master's  thesis  and  it  is  not  essential  for  the  first 
series  of  anticipated  tests. 

In  conjunction  with  these  physical  modifications  a  parallel 
development  of  the  LDA  data  acquisition  and  processing  circuitry  is 
in  final  stages  of  micropressor  development  involving  breadboarding 
etc.  Studies  directed  toward  processing  the  data  in  its  anticipated 
"sparse"  form  are  being  conducted  by  Dr.  Cole  and  his  graduate  (and 
undergraduate)  students.  This  work  is  described  and  recorded  in 
greater  detail  in  Appendix  A  to  this  report. 

Modelling  studies  of  a  preliminary  nature  have  been  completed 
for  the  particular  high  loading/unsteady  flow  phenomenon  described 
as  "propagating  stall".  A  very  simple  case  of  unstaggered,  uncambered 
thin  airfoils  has  been  analyzed  for  nonuniform  distribution  of  loading 
from  blade  to  blade  in  a  two-dimensional  cascade.  This  work  is 
detailed  in  Appendix  B  to  this  report.  The  principal  value  of  this 
latter  study  are  the  demonstration  of  the  possibility  of  lift-distribu¬ 
tion  anomalies  and  as  a  clear  guide  to  the  inclusion  in  the  model  of 
downstream  shed  vorticity  and  stability  considerations  which  do  not 
appear  in  the  work  shown  in  that  appendix. 

Personnel  who  have  been  supported  under  the  terminated  contract 
include,  in  addition  to  the  Project  Manager  and  the  Principal  Investiga¬ 
tor,  the  following  graduate  (G)  and  undergraduate  (U)  students  (the 
latter  being  assigned  exclusively  to  apparatus  construction  and 
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computing) : 


1 


Yao  Andrew  Tang  (G) 
Richard  Feltman  (G) 
Kenneth  Bird  (U) 
Stephen  Olson  (U) 


David  Sharp  (U) 
Raymond  Johnson  (U) 
Tom  Hughes  (G) 

Fady  Khairallah  (U) 


Stemming  from  the  nature  of  the  work  (design  and  procurement  of  an 
experimental  apparatus)  it  was  not  appropriate  to  assign  academic 
thesis  credit  to  the  graduate  students  for  their  contributions; 
it  is  anticipated  that  future  work  on  this  subject  will  involve 
both  analytical  and  experimental  phases  and  will  be  suitable  masters 
thesis  and  doctoral  dissertation  material.  The  elucidation  of  unsteady 
flow  phenomena  at  high  blade  loading  remains  a  fundamental  topic  of 
great  importance  in  designing  and  operating  axial  flow  compressors. 


CONCLUSIONS 

With  respect  to  data-reduction  methods  for  the  subject  LDA 
measurements,  several  conclusions  can  be  drawn.  From  an  extensive 
literature  search  and  review,  along  with  professional  consultations, 
it  is  concluded  that  the  theoretical  extensions  by  Papoulis  and  by 
Yen  are  apparently  the  only  existing  basis  for  dealing  with  the 
unusually  "sparse"  data  of  the  gated-LDA  situation.  Through  these 
extensions  or  adaptations  of  them,  it  is  theoretically  possible  to 
determine  unsteady  velocities  via  sparse  bursts  of  samples  randomly 
spaced  within  bursts.  From  numerical  simulations,  Papoulis'  "bunched- 
sample"  approach  is  concluded  to  be  satisfactorily  insensitive  to 
most  parametric  influences  on  the  data-reduction  process  (random 
error,  time  spacing  of  samples,  number  of  samples,  etc.).  However, 
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the  possibility  of  high  sensitivity  to  missing  data  samples  warrants 
further  study  before  this  approach  can  be  reasonably  assured  of 
viability.  Preliminary  investigation  of  Yen's  "migrated-sample" 
approach  however,  suggests  that  this  approach  may  not  be  satisfactory 
though  not  conclusively.  It  is  concluded  that  a  third  alternative 
justifies  further  study  via  simulation  -  Yen's  "minimum-energy" 
approach.  It  is  finally  concluded  that,  at  least  until  one  or  more 
of  these  approaches  can  be  more  definitively  found  satisfactory,  it 
is  appropriate  to  make  LDA  measurements  in  successive  blade  passages, 
interpreting  such  values  using  reasonable,  if  not  provable,  a  priori 
assumptions  concerning  the  nature  of  unsteady  velocities  in  neighboring 
blade  passages. 

With  respect  to  analytical/computer  model 1 ing  of  stall  "cells" 
at  high  cascade  loading  conditions,  a  beginning  study  has  shown  that, 
for  a  one-blade  cell,  boundary  conditions  may  be  satisfied  on  all 
blades  (stalled  and  unstalled).  The  solution  exhibits  local  regions 
of  reverse  flow.  It  appears  to  be  perfectly  straightforward  to  extend 
the  method  to  multiblade  cells  in  cascades  with  both  camber  and  stagger 
it  is  expected  that  a  richer  variety  of  individual  blade  loadings, 
within  and  without  the  cell,  will  be  exhibited  by  this  solution. 

The  features  which  are  lacking,  and  which  should  be  incorporated 
in  the  model,  in  any  future  work,  are  a  recognition  of  time  dependence 
and  consequently  the  induction  effects  of  the  shed  vorticity  downstream 
of  the  cascade.  With  the  successful  inclusion  of  these  two  key 
effects  in  the  model  it  should  be  capable  of  predicting  the  unsteady 
(vector)  velocity  at  a  fixed  point  in  the  cascade  flow  and  hence  the 
quantities  to  be  measured  by  the  LDA. 


APPENDIX  A 


GATED  DIGITAL  SAMPLING  FOR  TURBOMACHINES 
I.  INTRODUCTION 

The  prospect  of  LDA  measurements  of  in-rotor  air  velocities 
in  a  test  turbocompressor  leads  inherently  to  concern  in  the  area 
of  digital  data  processing.  This  Appendix  reports  results  of  a 
search  for  potentially  useful  theoretical  treatments  and  also 
subsequent  numerical  simulations  based  on  such  theory.  Simulations 
are  particularly  necessary  because  of  the  unconventional  situation 
presented  by  LDA  measurements  of  in-rotor  velocities.  The  desire 
to  measure  unsteady  in-rotor  velocities  presents  a  data-processing 
situation  apparently  undeveloped  heretofore. 

A  computer-aided  literature  search  uncovered  useful  existing 
theories  and  applications  of  pertinent  digital  signal  analysis. 

Because  of  the  inherent  randomness  and  the  sparseness  of  data 
sampling  which  was  expected  from  the  LDA/turbo-compressor  measurements 
particular  emphasis  was  given  in  the  literature  search  to  irregular 
or  non-uniform  spacing  of  samples  in  time.  Engineering  Index  and 
National  Technical  Information  Service  data  bases  were  searched 
resulting  in  a  variety  of  useful  references.  Over  30  references 
were  reviewed  as  a  result,  and  based  on  them,  further  secondary 
references  were  uncovered.  In  addition,  electrical  engineering 
and  mathematics  faculty  who  are  experienced  in  digital  sampling 
matters  were  consulted.  This  and  the  literature  search  led  to 
the  conclusion  that  the  subject  situation  was  not  dealt  with  effective 
ly  by  existing  developments.  A  review  of  the  conditions  of  conven¬ 
tional  digital  sampling  theory  makes  clear  the  distinctive  features 
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encountered  in  the  subject  LDA  application. 


It  is  well-established  in  the  theory  of  digital  signal 
analysis  ("sampling  theorem",  Whittaker,  1915;  Shannon,  1949) 
that  an  infinite  set  of  instantaneously-taken  samples  from  such 
a  signal  contains  all  the  information,  i.e.,  complete  specification, 
of  the  original  signal,  f(t),  iif  the  signal  is  bandwidth-limited, 
i.e.,  if  its  frequency  components  do  not  exceed  a  limiting  value 
(Stearns,  1975) .  In  fact,  the  signal  can  be  exactly  reconstructed 
from  its  sampled  values  by  appropriately  weighting  contributions 
from  all  the  sample  values: 


oo 

f(t)  =  L 

n=-oo 


_  .  sin[7rfs(t-tn) 

n  7Tfs(t_tn> 


( A- 1 ) 


where:  tn  is  the  time  of  the  nth  sample 

fs  is  the  frequency  of  sampling 

This  is  strictly  true  if  f(t)  is  limited  to  frequency  components 
less  than  1/2  fs  but  only  approximately  true  otherwise.  For  a  given 
bandwidth  -  limited  signal,  then,  the  sampling  frequency  must  exceed 
twice  the  highest  signal  frequency.  This  minimum  sampling  frequency 
is  sometimes  termed  the  "Nyquist  frequency". 

From  one  point  of  view,  the  need  for  limiting  the  bandwidth 
of  f(t)  can  be  thought  of  as  a  necessity  in  order  to  avoid  "aliasing", 
i.e.,  the  erroneous  interpretation  of  high-frequency  periodicity 
as  if  it  derived  from  signal  components  at  low  frequency.  Such 
incorrect  interpretation  is  familiar  hazard  of  using  a  stroboscope 
to  determine  rotational  speed. 
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A  second  source  of  error  in  the  practical  use  of  the  sampling 
theorem  arises  from  the  need  to  deal  with  a  finite  set  of  samples. 

In  practice,  "truncation"  of  the  theoretically  infinite  series 
representing  f(t)  leads  to  so-called  "truncation  error"  which  is 
unavoidable . 

Unfortunately,  however,  velocities  measured  by  an  LDA 
cannot  be  thought  of  as  even  a  finite  number  of  the  equally-spaced 
samples  deal  with  in  conventional  sampling  theory. 

II.  DIGITAL  SAMPLING  IN  LDA  MEASUREMENTS 

The  concept  of  laser  doppler  anemometry  has  several  alternative 
implementations,  the  most  common  method  can  be  visualized  as  small 
particles  (traveling  with  a  fluid  flow)  traversing  a  "measuring 
volume"  in  space  in  which  crossed  laser  beams  produce  a  constructive/ 
destructive  optical  interference  pattern.  As  particles  pass  from 
light  to  dark  regions  of  the  interference  "fringe"  pattern  within 
the  measuring  volume,  more  or  less  light  is  scattered  and  eventually 
received  by  a  photodetector.  This  results  in  a  time-wise  modulation 
of  the  photodetector  output  signal  related  to  the  speed  of  the 
scattering  particle.  The  time  at  which  such  a  particle  enters  the 
measuring  volume  is  essentially  random.  Therefore,  the  velocity 
value  determinable  from  the  modulated  photodetector  signal  must  be 
considered  as  representative  of  a  randomly  selected  time.  Successive 
determinations  are,  therefore,  not  at  the  equal  time  spacing  dealt  wi 
by  the  sampling  theorem  for  digital  signals  but  rather  at  random 
time  spacing. 

Despite  the  inherently-random  times  at  which  such  LDA  measure¬ 
ments  are  made,  it  is  conceivable  to  establish  corresponding 


equally-spaced  sample  values  to  which  the  sampling  theorem  and 
corresponding  data-reduction  methods  (e.g..  Fast  Fourier  Transform) 
could  apply.  This  might  be  accomplished  by  fitting  appropriate 
smoothed  curves  to  the  randomly  timed  data  samples  (as  an  approxi¬ 
mation  to  the  original  continuous  function  representing  time-varying 
velocities) .  These  smoothed  approximations  might  then  be  resampled 
at  equally  spaced  intervals  to  yield  the  requisite  sets  of  equally- 
spaced  sample  values  which  admit  conventional  data  reduction. 

Beside  resampling  in  order  to  give  equally-spaced  samples, 
it  is  even  theoretically  possible  to  deal  with  original  randomly- 
spaced  sample  values  themselves;  the  sampling  theorem  can  be  extended 
to  randomly-spaced  samples,  and  the  sufficiency  of  an  (infinite) 
number  of  such  samples  for  determining  signal  spectra  or  even 
reconstructing  signal  functions  is  assured  (Loo,  1969).  Thus,  if 
randomness  of  sample  spacing  were  the  only  complication  of  the  subject 
LDA/turbocompressor  measurements,  random-sampling  extensions  of 
conventional  theory  would  be  available  and  applicable.  Unfortunately, 
in-rotor  LDA  measurements  require  "gating"  of  the  usual  (randomly- 
spaced)  LDA  measurements,  i.e.,  accepting  only  those  measured  values 
occuring  within  well-defined  time  periods  as  is  discussed  in  the 
next  section. 

III.  GATED  DIGITAL  SAMPLING 

Measurements  of  gas  velocities  in  turbomachine  rotors  have 
been  accomplished  by  various  workers  who  have  determined  mean  and 
sometimes  root-mean- square  velocities  from  "gated"  LDA  measurements. 
This  approach  requires  the  accepting  of  data  only  if  it  is  measured 
during  a  brief  time  period  (gate  "open")  when  one  small  region 
within  one  blade  passage  falls  within  the  measuring  volume  of  the  LDA. 
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Figure  A-l  illustrates  this  scheme.  Successive  measurements  taken 
from  successive  revolutions  of  the  rotor  have  been  time-averaged  by 
some  researcners  and  r-m-s  values  calculated  in  order  to  "map" 
mean  gas  velocities  and  their  variances  or  deviations  within  rotor 
blade  passages.  To  move  beyond  mean  and  r-m-s  velocity  determinations 
into  unsteady,  time-resolved  measurements  forces  confrontation  with 
the  sampling  theorem  and  its  ramifications. 

If  only  (as  shown  in  Fig.  A-l)  a  single  velocity  value  is 
measured  during  each  gating  period  (once  per  rotor  revolution) , 
the  sampling  theorem  clearly  restricts  the  obtainable  information  to 
signal  components  of  frequency  less  half  the  rotor  speed  (half  the 
Nyquist  frequency) .  Components  at  such  low  frequencies  are  not 
of  greatest  interest  aerodynamically .  The  question  immediately  arises 
as  to  whether  it  is  possible  from  gated  LDA  measurements  to  provide 
information  concerning  higher  frequency  components  of  the  unsteady 
gas  velocities. 

Fortunately,  an  old  "folk  theorem"  of  digital  sampling  theory 
was  uncovered.  This  states  that  a  "signal  f(t)  may  be  represented 
by  anylinear  combination  of  irregularly  spaced  samples  f(tn),  provided 
the  average  sampling  rate  exceeds  the  Nyquist  rate,  i.e.,  that  the 
number  of  samples  per  unit  of  time  exceed  (on  the  average)  twice  the 
highest  frequency  present  in  the  signal  (Beutler,  1966).  Beutler 
has,  in  fact,  shown  that  "only  the  past  need  be  sampled  at  an  average 
rate  greater  than  the  Nyquist  rate  to  assure  error-free  recovery". 

In  this  context,  multiple  LDA  measurements  during  short  once-per- 
revolution  gating  periods  can  theoretically  provide  the  means  for 
deriving  measurements  of  unsteady  velocity  components  at  frequencies 
even  exceeding  the  rotation  il  freque  vcy  of  a  turbomachine  rotor. 
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However,  the  mathematical  theory  assuring  the  possibility  of 
LDA  measurements  of  unsteady  velocity  is  separate  from  the  question 
of  how  to  do  this  i.e.,  a  practical  data-reduction  scheme. 
Unfortunately,  the  literature  seems  to  be  devoid  of  a  data-reduction 
algorithm  or  approach  matched  to  the  gated  LDA  situation.  This 
forced  investigation  of  approacnes  available  for  similar,  if  not 
identical,  data-reduction  tasks  as  described  in  the  next  section. 

IV.  AVAILABLE  THEORY 

Two  available  approaches  to  digital  sampling  provide 
extensions  to  conventional  theory  which  are  distinctly  applicable 
to  the  subject  gated-LDA  measurements.  Other  authors  have  established 
the  sufficiency  of  samples  spaced  randomly  in  time  (e.g.,  Loo,  1969) 
and  the  significance  of  the  average  sampling  rate  as  the  equivalent 
(for  non-uniform  sample  spacing)  of  the  Nyquist-limit  sampling  rate 
(for  equally-spaced  samples)  (Beutler,  1966)  .  However,  it  appears 
that  only  the  two  works  cited  following  provide  expressions  for 
actually  carrying  out  the  theoretically- justified  reconstruction 
base'5  on  non-uniformly-spaced  sample  data. 

First,  Papoulis  (1977)  has  generalized  the  specifications  for 
sample  sets  from  which  a  signal  can  be  reconstructed  to  allow  for 
the  possibility  of  non-uniform  sample  spacing.  Papoulis'  formulation 
provides  for  "bunching"  of  samples,  i.e.,  a  non-uniform  sample 
spacing  within  a  finite-sized  "bunch"  of  samples  with  unending 
recurrence  of  the  same  non-uniform,  in-bunch  spacing  over  successive 
bunches  (see  Fig.  A-2) .  Reconstruction  of  the  original  signal,  f(t), 
from  sample  values  is,  in  general: 


oo 


f(t)  =  Z  gi (nT) y1 { t-nT)  +  g2 (nT) y2 ( t-nT) t  ...  +  g  (nT)y  (t-nT) 

n  =  —  oo 

(A-2) 

where:  T  is  the  period  of  recurrence  of  bunches 

n  is  the  running  index  designating  a  specific  bunch 
(of  non-uniformly-spaced  samples) 
m  is  the  running  index  designating  a  specific  sample 

(within  a  single  bunch  and  recurring  within  successive 
bunches ) 

g^  is  a  function  defining  the  samples  values  as  a  function 
of  time 

ym  is  the  reconstruction  function  for  all  samples  at  position 
"m"  in  all  bunches. 

Among  other  possibilities,  Papoulis'  theory  includes  pro¬ 
vision  for  sample  functions  g^(t)  which  take  the  value  of  the  signal 
at  some  well-defined  time  leading  or  lagging  the  sampling  time 
(i.e.,  for  arbitrarily  time-dispersed  sampling).  Alternatively, 
the  theory  allows  for  sample  functions  which  take  the  value  of 
one  or  more  derivatives  of  the  signal  function  as  well  as  the  function 
itself  at  the  sampling  times. 

Papoulis'  formulation,  then,  allows  generalization  of  the 
timing  and  definition  of  sample  values  but  always  within  the  context 
of  recurring,  repetitive  sample  bunches.  Figure  A-3  illustrates  two 
examples  of  this  generalization  for  which  Papoulis  developed  specific 
reconstruction  formulae.  The  second  example,  which  involves  bunching 
of  samples  in  twos,  can  be  seen  to  resemble  gated-LDA  sampling  but 
with  repeated  regular  sampling  patterns  rather  than  random  patterns 
during  successive  gates. 
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A  second  potentially  useful  extension  is  that  formulated  by 
Yen  (1956).  Yen  developed  four  specific  extensions,  one  of  which 
allows  for  the  "migration"  of  a  finite  number  of  conventional , 
equally-spaced  samples  to  arbitrary  sampling  times  (see  Fig.  A-4) . 

For  equally-spaced  samples,  reconstruction  of  limited- 
bandwidth  signal  functions  f(t)  from  an  infinite  sample  set  {f(tn)} 
is  accomplished  (as  described  above)  via  (after  Whittaker  and  Shannon) : 


f(t)  =  Z 

n=-oo 


sin  [ 7T  ( t-t  ) /T  ] 
f(tn)  TT(t-tn)/Tg 


(A-  3) 


where:  T  is  the  (fixed)  time  period  between  successive  samples, 

s 

Equally-spaced  samples  are  seen  to  result  in  the  same  weighting  or 

reconstruction  function  (sin  (ir  (t-t  ) /T  )/u(t-t  )/T  )  for  each 

ns  ns 

individual  sample  value,  f(tn). 

As  an  extension  for  non-uniformly-spaced  samples,  Yen 
considers  a  different  weighting  or  reconstruction  function  for  each 
different  non-uniformly  spaced  sample,  f(t  )  (total  of  N  values) 

r 

and  for  each  different  time,  t: 


N 

f(t)  =  Z  f(t)  Y  ( t)  +  Z  f(t  )  4»  (t) 

p=l  v  p  q — — oo  4  4 


non-uniformly-spaced  remaining  uniformly-spaced 

samples  samples 


(A-4) 


Yen  derives  general  expressions  for  the  (different)  4^  and  4^  functions. 

Another  of  Yen's  extensions  deals  with  a  finite  number  of 
arbitrarily-spaced  samples  used  to  reconstruct  the  original  signal. 

This  allows  unique  solution  if  the  resulting,  reconstructed  sample 
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function  is  constrained  to  minimize  its  "energy"  (integral  of  its 
amplitude  squared) . 

Both  of  Yen's  and  Papoulis'  extensions  are  significantly 
close  to  the  subject  LDA  situation  to  allow  dealing  with  gated  LDA  . 
measurements.  Papoulis'  recurrent  spacing  pattern  could  be  imposed 
by  smoothing  the  original,  randomly  spaced  LDA  sample  values  from  the 
LDA  and  resampling  (interpolating)  to  gain  the  necessary  recurrent 
pattern  of  sampling  times. 

However,  both  Papoulis*  and  Yen's  theories  apply  rigorously 
only  for  an  infinite  number  of  samples  (just  as  does  the  Whittaker 
reconstruction  algorithm  for  equally-spaced  samples) .  Neither 
addresses  the  question  of  errors  introduced  by  the  truncation  to 
finite-size  data  sets  as  required  in  practice.  A  further  necessary 
concern  for  error  which  is  lacking  in  Papoulis'  and  Yen's  treatments 
relates  to  the  specific  sort  of  bunching  of  data  anticipated  in  the 
gated  LDA  sutiation  -  sparse  data  occurring  in  bursts  with  relatively 
long  data-free  interludes  between  bursts.  No  treatment  of  such 
practically  crucial  concerns  was  found  in  the  literature.  Therefore, 
computer  simulations  were  carried  out  to  provide  insight  into  the 
practicality  of  applying  Papoulis'  or  Yen's  extensions  to  the  subject 
LDA  experiments.  These  are  described  in  the  following  section. 

A.  RECURRENT  BUNCHED  SAMPLES 

The  sampling- theory  generalization  by  Papoulis  (1977)  is 
particularly  easily  formulated  for  the  specific  case  of  two  data 
samples  per  bunch.  This  case  is  given  as  an  example  by  Papoulis. 

While  LDA  experiments  give  non-recurrent,  random  data  spacing  within 
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bunches,  such  data  could  be  converted  to  samples  satisfying  Papoulis' 
example- scheme  by  data  smoothing  and  resampling.  A  computer  simulation 
was  programmed  in  order  to  evaluate  parametric  influences  on  the 
accuracy  of  signal  reconstruction  and  subsequent  spectral  analysis 
according  to  Papoulis'  scheme. 


A. 1  INFLUENCE  OF  RELATIVE  FREQUENCY  OF  SIGNAL  &  SAMPLING 


Several  simulation  runs  were  made  with  sample  sets  of  6  pairs 
(bunches)  of  samples  taken  from  a  pure  sine  wave.  To  approximate 
the  sparse  data  of  the  grated  LDA  situation  the  time  interval 
between  the  paired  sample  of  each  "bunch"  was  0.02  units  on  a  scale 
of  1  time  unit  between  bunches  (see  Fig.  A-5) .  Figure  A-6  shows 
the  original  (sinusoidal)  signal  and  the  pairs  of  sample  points 
superimposed  for  a  case  in  which  the  signal  frequency  is  0.9  of  the 
frequency  at  which  bunches  of  samples  are  taken.  The  frequency 
of  occurrence  of  bunched  samples  corresponds  with  the  gating 
frequency  for  gated-LDA  measurement  (or  the  turbomachine  rotational 
speed  in  the  subject  experiments) .  For  the  purpose  of  these  and 
further  plots,  times  are  non-dimensionalized  by  the  period  between 
gates  (period  of  turbomachine  rotation) . 


Figure  A-7  shows  a  reconstruction  of  the  signal  of  Figure  A-6. 
The  reconstruction  is  based  on  the  (bunched)  sample  points  (which 
are  shown)  and  on  Papoulis'  reconstruction  equation: 


2  7T  2tt  , 

COSjjT-a  -  COSjjT-t 


f(t)  = 


2ir 

T 


.  2u 
sin^r-a 


CO 

E 

n  =  -a> 


f (nT+a) 
t-nT-a 


f (nT-a) 
t-nT+a 


] 


(A-5) 
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where  f(t)  is  the  signal  function 


2a  is  the  (non-dimensional)  time  interval  between  paired 
samples  (within  a  burst) 

T  is  the  (nominal)  sampling  period  (time  between  sample 
bunches) 

n  is  the  running  index  designating  different  bunches 

The  figure  illustrates  the  occurrence  of  "leakage",  i.e., 
non-zero  values  (oscillations)  of  the  reconstructed  signal  outside 
the  time  span  of  the  original  sampling  (approx,  t  =  0  to  t  =  5)  . 

The  figure  also  shows  that  reasonably  faithful  reconstruction  is 
feasible  even  when  the  signal  frequency  (0.9  in  Figs.  A-6  and  A-7)  is 
close  to  the  nominal  bunched-sample  frequency  (1.0).  Similar  results 
are  shown  in  Figure  A-8  for  which  the  relative  phase  of  the  signal 
and  the  bunched  sampling  is  varied,  maintaining  the  same  frequencies 
as  in  Figures  A-6  and  A-7. 

Successively  increasing  the  signal  frequency  to  well  above 
the  nominal  sampling  frequency  shows,  first,  increased  distortion 
and  finally  gross  error  in  the  reconstructed  signal.  This  trend 
reflects  the  theoretical  limitation  of  the  reconstruction  process 
(even  using  an  infinite  sample  set)  to  signal  frequencies  less  than 
one-half  the  time-averaged  sampling  frequency.  In  this  case,  the 
time-averaged  sampling  frequency  is  (2  samples  per  bunch)  x  (1  time 
unit  per  bunch)  =  2  samples  per  time  unit.  Therefore,  the  signal 
frequency  should  not  exceed  1.0  cycle  per  time  unit.  In  fact,  the 
reconstruction  of  signals  of  frequencies  1.1  and  1.3  cycles  per  unit 
time  (Fig.  A-9)  is  seriously  in  error. 
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INFLUENCE  OF  SAMPLE  MEASUREMENT  ERRORS 


The  effect  of  random  experimental  errors  in  measuring  sample 
values  was  also  simulated.  Random  errors  were  added  (or  subtracted) 
to  the  calculated  (sinusoidal-signal)  sample  values  which  were 
previously  used.  The  random  errors  were  taken  to  have  zero  mean 
value  and  an  arbitrarily  specified  standard  deviation.  Fig.  A-10 
illustrates  the  increasing  degradation  of  the  reconstructed  signal 
in  terms  of  reconstructed  signal  shape  as  the  standard  deviation  of 
experimental  error  increases.  At  least  for  the  baseline  conditions 
of  this  simulation  (12  sample  points,  signal  freq.  =  0.9  max.  allow¬ 
able  signal  freq..  Figs.  A-6  and  A-7)  ,  random  error  with  standard 
deviation  exceeding  several  percent  is  clearly  unacceptable.  Even 
doubling  the  number  of  samples  does  not  seem  to  alter  such  a 
conclusion  much  (Fig.  A-ll)  as  is  especially  evident  in  the  resulting 
power  spectra  (Fig.  A-12) . 

A. 3  INFLUENCE  OF  NUMBER  OF  SAMPLES  (SAMPLING  DURATION 

Simulations  for  varying  numbers  of  bunches  of  paired  samples 
were  carried  out.  Fig.  A-13  illustrates  the  results,  with  respect 
to  shape  of  the  reconstructed  signal  itself  and  Fig.  A-14  shows  the 
corresponding  power  spectra.  Increased  numbers  of  data  are  seen 
to  improve  spectral  discrimination  considerably  as  is  expected,  but 
the  high  spectral  signal/background-noise  (e.g.,  7/1,  Fig.  A-14) 
for  a  modest  number  of  samples  (24)  was  unexpectedly  encouraging. 

A. 4  INFLUENCE  OF  GATE  DURATION 

In  the  subject,  gated-LDA  situation,  the  "gate-open"  interval 
is  but  a  small  fraction  of  the  total  experiment  (sampling)  interval, 
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for  example,  0.001  to  0.01.  Short  "gate-open"  intervals  imply 
decreased  probability  of  a  valid  LDA  measurement  during  that 
interval  but  increased  spatial  resolution  in  measuring  turbomachine 
rotor  blade-passage  velocities.  Short  gate-open  intervals  can  also 
be  expected  to  degrade  the  accuracy  of  reconstructing  signals  or 
gaining  spectral  data  from  samples  taken  during  such  short  intervals. 

To  help  evaluate  these  effects,  several  simulations  computer  runs 
were  made  with  varying  time  between  the  (paired)  samples  within  a 
bunch  or  burst  (2a  in  Eq.  A-5) . 

Fig.  A-15  and  A-16  show  how  reconstructed-signal  shapes 
and  their  spectra  respond  to  increasing  bunching  of  samples  without 
change  in  average  sampling  rate.  Clearly,  very  sparse  samples,  (i.e., 
widely  separated  bunches  or  "bursts")  are  feasible  without  severe 
degradation  of  signal-reconstruction  possibilities  using  Papoulis' 
scheme. 

B.  FINITE  SET  OF  NON-UNIFORMLY-SPACED  SAMPLES  (EMBEDDED  IN  INFINITE 
SET  OF  UNIFORMLY-SPACED  SAMPLES) 

The  sampling-theory  generalizations  by  Yen  (1956)  appear 
particularly  applicable  to  the  subject  gated-LDA  measurements. 

Unlike  Papoulis'  generalization,  those  of  Yen  allow  arbitrarily- 
spaced  samples,  consistent  with  the  gated  LDA  measurements  which 
are  expected  to  occur  at  random  times  within  sparse  bursts. 

Yen's  approach  has  been  demonstrated,  on  the  other  hand,  to  require 
relatively  long  computation  times  (Sankur  and  Gerhardt,  1973). 

This  is  not,  however,  a  crucial  shortcoming  for  the  subject  application. 


Much  more  likely  to  be  practically  significant  is  the  fact  that 
Yen's  reconstruction  functions  can  be  expected  to  assume  large 
values  between  sample  points.  This  leads  to  potential  computational 
difficulties.  The  required  summation  of  terms  (Eq.  A-4)  seems 
necessarily  to  lead  to  calculating  small  differences  between  large 
numbers  (from  different  terms)  yielding  potentially  large  errors. 

This  potential  problem  as  well  as  the  problem  of  errors  due  to 
truncation  to  a  finite  number  of  samples  has  led  to  numerical  computer 
simulation.  Computer  programming  to  allow  this  simulation  is  in 
progress  but  at  this  date  has  not  yet  been  completed,  and  its  results 
cannot  be  reported  here.  Preliminary  results  do,  however,  suggest 
strongly  that  the  problems  cited  above  are  probably  excessive 
rendering  this  alternative  impractical. 


C.  NON- UNIFORMLY-SPACED  SAMPLES  AND  "MINIMUM  ENERGY"  RECONSTRUCTION 
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BUNCHES  OF  IRREGULARLY-SPACED  SAMPLES 
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FIGURE  A- 2;  Schematic  of  Papoulis'  "Bunched"  Samples 


(b)  Samples  of  Function  Values  in  "Bunches"  of  Two 
FIGURE  A-3:  Examples  of  Non-Conventional  Sampling  Schemes 
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FIGURE  A-4:  Schematic  of  Yen's  "Migrated"  Samples 


FIGURE  A-5:  Sampling  Scheme  for  Simulations 
(Not  to  Scale) 


FIGURE  A- 6 ;  Original  (Sinusoidal)  Signal  Function 
Showing  Bunched  (Paired)  Samples 
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FIGURE  A- 7 ;  Signal  Reconstructed  from  Bunched  Sample 


FIGURE  A-9  (CONTINUED) 


STD.  DEV.  OF  ERROR  =  0.01 


FIGURE 
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STD.  DEV.  OF  ERROR  =  0.05 

-10 :  Signal  Reconstructions  for  Various 

Experimental  Errors 


NON  Din. FREQUENCY 
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STD.  DEV.  OF  ERROR  =  0.05 


FIGURE  A-12:  Spectra  for  Reconstructed  Signals 
.  -  with  Error  and  Increased  Number 

of  Samples  (48) 
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a  =  0.10  (•),  a  =  0.001  (0,2  pts.) 

FIGURE  A-15;  Signal  Reconstructions  for  Various  \ 

Gate-Open  Intervals  ] 

(a  =  Non-Dim.  Half  Interval)  j 
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FIGURE  A-16;  Spectra  of  Reconstructed  Signals 
For  Various  Gate-Open  Intervals 

(a  =  Non-Dim.  Half  Interval) 


Appendix  B 

CASCADES  OF  NONUNIFORM  LOADING 

We  wish  to  study  initially  steady  nonuniform  distributions 
of  lift,  or  circulation,  along  cascades  of  airfoils.  The  simplest 
problem  of  this  nature  consists  of  an  unstaggered  infinite  cascade 
of  flat  plate  airfoils  in  two  dimensional,  incompressible  flow. 
Consider  the  case  where  alternate  blades  of  unit  semichord  have 
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the  same  distribution  of  vorticity, 
and  hence  only  two  families  of 
interdigitated  blades  may  be  dis¬ 
tinguished.  Referring  to  the  sketch, 
the  vertical  velocities  induced  on 
the  chordlines  by  this  system  may 
be  shown  to  be  given  by 


U 


II 


vI(x)  =  ~  K(£-x)YlU)d£;  +  J  L(5-x)YlI(C)d? 

vn(x)=  —  K(£-x)  YlI  (C)d?  +  Jg  f  L(C-x)  Y][  (C)dC 


(1) 


(2) 


where 


K(C-x)  =  ~  coth  Ijg  (C-x)] 


(3) 


L(C-x)  =  tanh  [^r  (£-x)]  (4) 

It  is  the  distributions  y  and  y.^  that  are  sought  in  this 
investigation.  The  kernels  K  and  L  are  correct  only  for  the  present 
problem  of  an  alternating  distribution  in  passing  from  blade  to  blade; 
their  derivation  is  perfectly  straightforward.  The  largest  magnitude 
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of  the  argument  of  the  hyperbolic  functions  that  can  be  expected 

is  about  3  or  4,  hence  the  kernels  can  be  approximated  quite  accurately 

by 

K  (x-£)  =  K-x)'1  +  Aj^S-x)  +  a2(£-x)3  (3*) 


L(x-C)  = 


B1(C-x)  +  B2(S-x) 


where  the  constants  A^,  A2,  B^  and  B2  may  be  fitted  by  least 
squares  over  the  interval  0<  £-x  < tt  (2s)  Note  that  the  pitch/chord 

ratio  is  s/2. 

To  solve  the  coupled  integral  equations  (1)  and  (2)  it  is 
convenient  to  transform  the  variables  according  to 

5  =  cos0  ,  x  =  cos*  ,  0  £  0  or  if  <  it  (5 

In  the  arguments  0  and  4>  the  distributed  circulations  may  be  expanded 
in  the  so-called  Glauert  series 


(0)  =  t-j.  cot  *50 

+  Xj  tan  *50 

+  I 
1 

oo 

Cln  sin 

n0 

(6) 

(0)  =  TjjCOt  *50 

+  X^tan  *50 

+  E 
1 

Cllnsin 

n0 

(7) 

These  distributions  produce  circulations  round  the  representative 
airfoils  equal  to 

rifII  =  71  ^  TI , II  +  XI,II  +  *  CI1,II1) 

With  these  substitutions  equations  (1)  and  (2)  may  be  transformed 
as  follows 

V.j.(*)  =  2-jj-^[(cOS0  -  cos*)  +  A^  (cos0  -  COS*)  +  A2(COS0  -  cos*)J] 
o 


[tt  +  X  +  (tt  -  X_)cos0  +  £  CT  sin0  sin  nO]dO 
ill!  ^  In 

TT 

+  yjj-  J  [B^(cos0  -  cos*)  +  B2(cos0  -  cos*)3]  . 
o 

00 

[(in  +  Xjj)  +  cos©  +  i  CIIn  sin0  sin  n0]dO(l') 


1  11  -1  3 

vI I  ( 4> >  =2 ~  £  [  (cosO  -  cosO)  +  (cosO  -  cos$)  +  A2(cos0  -  cos 4>)  ] 

o 

00 

tin  +  Xjj  +  (ijj  -  A  )cosO  +  E  CIIn  sinGsin  nO]dQ 

IT  ^ 

+  J  [B^  (cosG  -  costfi)  +  B2  (cosO  -  cos<J>)^] 
o 

OO 

[xT  +  A  +  (t_  -  A_)  cosG  +  E  CT  sinO  sin  nO]dO  (2*) 

I  j.  I  I  ^  In 

Using  the  well  known  singular  integral 

TT 

f  cos  nQ  _  sin  n<J>  , „ 4 

o  cosO-cosp  sit\(p 

the  various  trigonometric  integrals  may  be  evaluated  with  the 

following  result. 

2Vj.  ( <f>)  =  (Tj  -  Aj)  (1  +  2  A1  +  8  A2  ^  +  CI2*T  A1  +  2  A2*  +  CI4*l6  A2* 

+  (TII  "  XII) (2  B1  +  8  B2)  +  CII2(4  B1  +  2  B2)  +  CII4(l6  B2) 

+  cos<(>(-CI1  -  (tj.  +  *3;  +  2  CI1)A1  +  (~4  TI  "  1  XI  ~  4  CIl  “  8€I3)A2J 
+  cos4>[-(tii  +  AII  +  2  CII1)B1  +  (~4  TII  _  4XII  "  4  CII1~ICII3)B21 
+  cos2<|>[-Ci2  +  <4  T!  ~  4  +  g  Ci2)A2^  +  cos24>  I  (4TII  ~  '4XII+8CII2)B2? 

+  cos3<})[-Ci3  +  (4  Tj  +  4  +  g  Cil^A2^  +  cos3<j>  [  B2^ 


-  E  CIn  cos  n  <p  .  (10) 

4 

2vtt  (4* )  —  (  “  ^11  ^  ^  "o  o'  )  +  ^1  *-TTii  ( TT  A0 )  , 


2  18  2'  112 ' 4  1  2  2‘ 


114  16  2' 


etc . 


i.e.  the  subscripts  I  &  II  are  interchanged  throughout  2v^.  to  obtain  2v^j. 

The  boundary  condition  on  the  flat  plat  airfoils  requires 
that  the  induced  velocities  on  the  chordlines,  Vj  and  v^,  be  equal 
and  opposite  to  the  onset  velocity  component  V  of  the  mean  flow,  a 


constant.  The  expression  of  this  B.C.  leads  to  the  matrix  equation 
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where 

[A]  {C}  =  -2VU0001000} 

C  is  the  vector  of  coefficients 

{C}  =  {t z  Xr  C 

11  CI2  CI3  TII  XII  CI1 

and  the  8  x  1 
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In  terms  of  the  4x5  submatrices  [a]  and[b]  it  is  clear  that 


where  the  elements  of  [a]  and  {b]  may  be  determined  from  [A] 
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The  solution  of  these  equations  represents  the  possible 
configurations  or  distributions  of  circulation  (bound  vorticity) 
on  the  two  families  of  airfoils  constituting  the  cascade.  To 
perform  a  systematic  study  of  these  solutions  it  is  convenient 
to  consider  the  four*  cases  delineated  in  the  following  table 


Case 


II 


II 


1 

ii 

iii 

iv 


0 

0 


0 

0 

0 


Afterwards  a  more  general  study  may  possibly  be  attempted. 

Case  (i)  [t^  =  T  =  0].  This  is  the  normal  steady  flow  case, 

evidently  stable  in  time  and  with  the  expected  results:  A  1  =  A 
C-j- 1  =  Cjj^,  etc.  Both  families  of  airfoils  have  identical  distribu¬ 
tions  and  hence  the  cascade  loading  is  uniform  in  that  sense.  By 
making  use  of  this  demonstrable  conclusion  the  equations  may  be 
simplified  with  the  solution  easily  obtainable. 

Upon  deleting  the  1st  and  6th  columns  of  matrix  (14)  and  the 

1st  and  6th  rows  from  the  vector  (13),  it  is  found  that  equation 

(12)  is  separable  into  two  equations  that  may  be  solved  successively 

i.e.,  first  for  A  A  T  . 
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*Two  other  cases  may  be  discerned,  but  those  case 
the  roles  of  I  &  II.  n_c 
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0 


*■*■»'*  >»*  i**ct 


The  symmetry  displayed  in  the  coefficient  matrix  of  equation  (15) 
results  in  A^  =  A.^,  CI2  =  CII2  anc*  hence  the  equations  may  be  reduced 
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bn 

a13  +  b13 

31 

+  bn 

a33  +  b33 

=-2V 


where  A  =  A^  or  A^  and  C 2  =  C I2  or  CII2.  The  coefficient  matrix  in 
equation  (16)  may  also  be  partitioned  symmetrically  into  submatrices 
with  the  conclusion  that  =  C  CI3  =  CII3  and  hence 


22 

+  b22 

a24  +  b24 

42 

+  b42 

a44  +  b44 

J 

=  -A 


a21  +  b21 


a41  +  b41 


where  =  C.^  or  Ci;rl  and  =  CI3  or  CII3  . 

A  rougher  approximation  may  be  obtained  by  letting  A2  =  b2  =  0 
in  which  case  it  is  found  that  =  CII2  =  =  CII3  -  0  and 

A  +  B 

XI  "  XII  "  1+JjA1+J5B1  '  CI1  =  CII1  =  ~X  1+J5A1+%B1 


Finally,  to  recover  the  familiar  isolated  airfoil  result  it  is 
only  necessary  to  let  A.  =  B,  =  0  and  so  arrive  at 


til 


The  role  of  U,  a  measure  of  the  massflow  or  throughput. 


determines  the  incidence  for  specified  V,  and  hence  the  cascade  loading 
flow  deflection,  etc.  In  order  to  explore  the  local  flow  magnitudes, 
stagnation  points  and  possible  flow  reversals  it  is  interesting  to  take 
a  short  digression  and  consider  the  flow  pattern  associated  with  a 
column  of  point  vortices  of  strength  r  =  tt(A  + 

Column  of  Vortices 

Each  airfoil  has  a  circulation  whose  strength  may  be  expressed 
as  T  =  it  (A  +  =  2ttVK,  where  X  and  in  turn  are  proportional 

to  V  so  that 


K  =  ^  Mv)  +  5sC1(V) 
V 


(8i] 


A  vertical  row  of  stationary  vortices  of  strength  T  =  2ttVK 
and  spacing  s  immersed  in  an  onset  flow  with  components  U,  V  has  a 
complex  potential  (z  =  x  +  iy) 

F  =  <J>  +  ill/  =  (U  -  iV)z  +  iVK  log  sinh  ^  (19] 

The  corresponding  stream  function  is 

il>  =  Uy  -  V [x-!jK  log  (sinh2  ~  +  sin2  ^  ))  (20] 


and  the  complex  velocity  is 

u  -  iv  =  U  -  iV  +  — —  coth 

s  s 

Extracting  the  u,v  components  results  in 

sin  (27iy/s) _ 

cosh  ( 2ttx/s )  -  cos(2iTy/s) 

sinh  ( 27tx/s) _ 

cosh  (2tjx/s)  -  cos(2Tiy/s) 


u  =  U  +  -  VK 
s 

and 

v  =  v  -  Z  VK 
s 


(21) 
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Analysis  of  Eq.  (22)  shows  that  a  reverse  flow  region  may  exist 
on  the  intervals  n  -  <y/s  <  n  for  any  integer  n  and  for  x  near 
zero.  This  reversal  region  is  enlarged  for  small  values  of  U 
relative  to  V,  i  ,e.  for  large  incidence.  As  such  it  may  be 
thought  of  as  a  "stall"  zone,  although  it  occurs  near  the  nose 
on  the  pressure  side  of  the  airfoils  which  are  modelled  by  the 
row  of  discrete  lifting  vortices.  Consequently,  this  type  of 
reverse  flow  should  be  associated  with  "cascade  stall".  The 
nonoccurrence  of  such  reverse  flow  patches  may  be  explained  by 
stability  considerations  which  dictate  a  shift  to  another  flow 
pattern  as  the  incidence  is  progressively  increased. 

Column  of  Countervortices 

A  second  digression  is  instructive  at  this  point.  A  column 
of  equispaced  vortices  of  constant  strength  T  =  2ttVK  but  of  alter¬ 
nating  sign  is  the  simplest  analog  of  case  (ii)  to  follow. 


The  complex  potential  for  this  case  is 
F  =  4>  +  iip  =  (U  -  iV)z  +  iVK  log  tanh  ^ 
with  corresponding  stream  function 

U»  =  Uy  -  V  (x  -  |K  log  cosh  —  +  cos  } 

,  iTX  iry 

cosh  —  -  cos  — 
s  s 


(23) 


(24) 


The  u,v  components  of  velocity  for  this  flow  are 

u  =  U  +  ^ VK  sinny/s  cosh  ttx/s _ 

s  cosh^  ttx/s  -  cos^  Tiy/s 


and 


v  =  V  -  n—  cosTty/s  sinh  ttx/s _ 

s  cosh^rx/s  -  cos^  Tty/s 


(25) 


(26)  J 

I 


*It  is  instructive  to  trace  the  streamlines  ty=0,  <p=Us,  ijj=2Us,^=3Us. 
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This  gives,  near  x=0,  a  sequence  of  horizontal  "jets",  alternating 
between  right-  and  left-flowing  and  with  v-components  to  satisfy 
continuity.  The  velocity  U  biases  these  jets  to  have  larger  core 
velocities  in  the  right-flowing  "jets".  The  perturbations  die  out 
with  |xj  >  0. 


Case  (ii)  [t^O,  X^.^0]  •  This  requirement  of  a  Kutta  condition 
(smooth  flow)  at  the  leading  edge  of  family  -  II  airfoils  should 
provide  additional  insight  concerning  reversed  flows.  Displaying 
the  matrices  that  obtain  for  this  case 
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Comparing  equations  (17)  and  (15)  it  is  clear  that  the  sign  changes 
in  the  third  column  imply  in  case  (ii) 


(28) 


"TII  XI'  CI2  “  CII2 

and  equation  (15*)  (in  which  X  is  identified  with  X^)  may  be 
employed  to  find  these  values. 


(29) 


The  symmetry  conditions  applied  to  equation  (18)  result  in 


CI1  -CII1'  CI3  ~CII3 


(30) 
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and 


-  _  -1 

a22  ~  b22  a24  “  b24 

CI1 

a21  b21 

a42  “  b4 2  a44  "  b44 
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=  ~ 

a41  "  b4 1 

The  most  interesting  conclusion  of  studying  case  (ii)  is  the  fact 

that  Tj-j.  is  negative,  presuming  A  is  positive.  This  means  that  there 

is  flow  rearward  along  the  pressure  surface  of  the  family-II  profiles 
near  the  trailing  edge.  Questions  of  separation  at  the  sharp  trailing 

edge  are  ignored  for  the  time  being.  The  flow  then  turns  around 

the  TE  and  flows  forward  (upstream  in  the  sense  of  the  mean  flow) 
along  the  suction  surface.  Since  A^.  =  -Tjj,  =  -C  ^  and  = 

Tj  =0,  the  circulation  alternates  between  ±  tMA^  +  h  C.^)  from  blade 
to  blade;  there  is  no  net  deflection  of  the  throughflow  far  from 
the  cascade.  This  could  be  used  as  a  definition  of  "cascade  stall". 

Case  (iii)  [1^=0,  A^^O]  .  There  is  no  solution  because  C  is 
singular  for  this  case.  Physically  this  corresponds  to  the  fact 
that  "shock-free"*  flow  cannot  be  obtained  with  an  unstaggered 
cascade  of  uncambered  airfoils,  except  at  the  unloaded  (zero  deflection) 
condition.  Evidently  this  restriction  carries  over  to  nonuniform 
cascade  flow  when  only  alternate  blades  are  supposed  to  be  shock-free. 

Case  (iv)  [Aj  =  A  =0]  •  The  solution  for  this  case  may  be  deduced 
by  considering  a  reflection  of  case  (i) .  Thus  =  t  =  negative 
number  is  an  immediate  conclusion.  Comparison  with  equations  (15) 
and  (15 ')  leads  to  the  solution 


"  all  +  bll 

a13  +  b13 

m  m 

-T 

=  -2V 

’  1 

_  a31  +  b31 

a33  +  b33  _ 

_C2 

0 

*  local  incidence  equals  zero  at  the  leading  edge. 


and  similarly,  a  study  of  equations  (16)  and  (16  )  leads  to 
the  conclusion 
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where  t  =  ij  or  TII'  Cl  =  CH  or  Cnl'  etc* 


(33] 


Since  the  circulation  around  each  airfoil  is  now  (case  (iv)) 

r  =  1T(T  +  h  Cj)  (34 

and  this  is  a  negative  number,  the  airfoils  of  case  (iv)  produce 
a  negative  deflection  of  the  flow;  the  mean  flow  is  accelerated  in 
passing  through  the  cascade. 

Conclusions 

Cases  (i) (ii)  and  (iv)  taken  together  present  three  regimes 
of  cascade  flow  having  positive,  zero  and  "negative"  flow  deflection. 
When  the  horizontal  component  of  mean  velocity  U  is  considered. 

Case  (i)  corresponds  to  normal  flow  from  left  to  right  at  low  incidence; 
case  (ii)  is  a  kind  of  "cascade  stall"  in  which  the  deflection  drops 
to  zero  and  case  (iv)  is  best  interpreted  as  reverse  flow  with 
uniform  loading  along  the  cascade.  Plotting  upstream  and  downstream 
flow  angles  &  02  as  a  function  of  V/U  as  in  the  sketch,  illustrates 
this  behavior  (unit  solidity  assumed) . 


1 


The  transition  from  type  (i)  to  type  (ii)  characteristics  is 
unknown.  In  fact  stability  considerations  must  be  invoked  with 
the  likely  result  that  case  (ii)  is  not  an  attainable  steady  state 
configuration.  A  purely  hypothetical  transition  is  shown  by  the 
solid  lines  in  the  sketch. 

Stability 

In  order  to  study  stability  one  must  allow  for  time  dependency 
by  allowing  the  X,  x  and  C  coefficients  to  vary  with  time,  and  also 
accounting  for  the  shed  vorticity  as  the  circulations  Tj,  change 

with  the  evolution  of  time.  This  analysis  will  be  attempted  and 
reported  in  a  future  report. 

Example 

A  numerical  example  was  pursued  for  the  following  parameter 
values: 
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